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ABSTRACT 

We present Spitzer infrared observations to constrain disc and dust evolution in young T Tauri stars in MBM 12, a star-forming cloud 
at high latitude with an age of 2 Myr and a distance of 275 pc. The region contains 12 T Tauri systems, with primary spectral types 
between K3 and M6; 5 are weak-line and the rest classical T Tauri stars. We first use MIPS and literature photometry to compile 
spectral energy distributions for each of the 12 members in MBM 12, and derive their IR excesses. Of the 8 stars that are detected 
with MIPS (spectral types between K3 and M5), only 1 lacks an IR excess - the other 7 all have an IR excess that can be attributed 
to a disc. This means that in MBM 12, for the detected spectral types K3-M5, we have a very high disc fraction rate, about 90%. 
Furthermore, 3 of those 7 excess sources are candidate transitional discs. The four lowest-mass systems in the cloud, with spectral 
types of M5-M6, were undetected by Spitzer. Their upper limits indicate that they either have a transitional disc, or no disc at all. The 
IRS spectra are analysed with the newly developed two-layer temperature distribution (TLTD) spectral decomposition method. For 
the 7 T Tauri stars with a detected IR excess, we analyse their solid-state features to derive dust properties such as mass-averaged grain 
size, composition and crystallinity. The mass-averaged grain size we determine from the 10 micron feature has a wide range, between 
0.4 and 6 pm. This grain size is much smaller in the longer-wavelength region: between 0. 1 and 1.5 pm. We find that later-type objects 
have larger grain sizes, as was already shown by earlier studies. Furthermore, we find a wide range in mass fraction of the crystalline 
grains, between 3 and (at least) 30%, with no relation to the spectral type nor grain size. We do find a spatial gradient in the forsterite 
to enstatite range, with more enstatite present in the warmer regions. The fact that we see a radial dependence of the dust properties 
indicates that radial mixing is not very efficient in the discs of these young T Tauri stars. The sources that have the least amount of 
disc flaring have the largest grain sizes, pointing to dust settling. A comparison between the objects with companions closer than 400 
AU ('binaries') and those with wider or no companions ('singles'), shows that disc evolution already starts to differentiate between 
both groups at an age of 2 Myr: the excess at 30 pm is a factor 3 larger for the 'single' group. The SED analysis shows that the discs 
in MBM 12, in general, undergo rapid inner disc clearing, while the binary sources have faster disc evolution. The dust grains seem 
to evolve independently from the stellar properties, but are mildly related to disc properties such as flaring and accretion rates. 

Key words. Stars: pre-main sequence; (Stars:) planetary systems: protoplanetary disks; Infrared: stars; Stars: late-type 



1. Introduction 

Our understanding of the star and planet formation process has 
advanced through comprehensive observations of nearby star- 
forming regions, such as the Taurus-Auriga cloud or the much 
denser and massive Orion Nebula Cluster. Because these sites 
are relatively close and young, it is often possible to determine 
the complete census of cluster members, from the highest mass 
stars down to brown dwarfs with masses of just a few Jupiters 
(e.g. Luhman et al. 2003 Muench et al. 2002 ). Furthermore, the 
detection of companions through adaptive optics and other high- 
resolution imaging and the subsequent derivation of multiplicity 
rates is an important piece of the star formation (SF) puzzle (e.g. 
Petr et al. I1998I Ratzka et al. |2005l >. In addition, near-IR stud- 
ies of SF regions have shown that young objects exhibit near-IR 
excesses, commonly attributed to the presence of protoplanetary 
discs (e.g. Adams et al. 1987, Lada et al. 2000). The structure, 



mass and size of those discs provide constraints on the process 
of planet formation in these young environments. 

To study the young cluster members in more detail, it is im- 
portant to characterise their circumstellar environment, which is 
best done at infrared and/or millimetre wavelengths. With the 
launch of Spitzer, this field has made huge steps forward, as for 
the first time, it was possible to observe the disc and dust char- 
acteristics of large samples of T Tauri stars (e.g. Kessler-Silacci 
et al. 120061 Sicilia-Aguilar et a l. 120071 Pascucci et al. 120081 ) and 
even brown dwarfs (Apai et al. 2005 , Scholz et al. 120071 ). 

Cloud 12 in the catalogue of Magnani, Blitz and Mundy 
(119851 ), in short MBM 12, is a molecular cloud at high galac- 
tic latitude (/, b - 159.4°, -34.3°), with relatively high extinction 
(Ay > 5 mag). Radio maps of the region in CO give a mass esti- 
mate for the whole complex of 30 to 200 M Q (Pound et al. |19901 >. 
The cloud does not appear to be gravitationally bound, and the 
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molecular gas is expected to dissipate within the next few Myrs 
(Zimmermann & Ungerechts I1990I I. Hearty et al. (2000) used 
ROSAT observations to detect X-ray sources in MBM 12, and 
followed-up stellar candidates by optical spectroscopy, to con- 
firm 8 X-ray emitting T Tauri stars with an upper age of 10 Myr. 
It is one of the few clouds at high galactic latitude known to 
harbour such young objects. From its content, the cloud can be 
seen as a precursor of a TW Hydrae-like association, where the 
molecular material has not yet disappeared. Although MBM 12 
was initially thought to be one of the nearest star formation re- 
gions, at a distance of only 65 pc, it is now determined to lie at a 
distance of 275 pc (Luhman l2001l l. 

A census of this association - complete down to masses of 
0.03 M - was determined by Luhman (2001), based on sen- 
sitive near-IR and optical photometry. The candidate members 
were confirmed by follow-up spectroscopy, which also allowed 
to study their Ha and Li 6707 A line properties. In total, 12 T 
Tauri stars were found to be real members of MBM 12, and an 
age of 2^ Myr was derived from both the lithium line and the 
location of the objects in the H-R diagram. The spectral types 
for those 12 TTS is between K3 and M6. An additional study 
by Luhman & Steeghs (2004) of 7 candidate members could not 
confirm any more members. The T Tauri stars in MBM 12 have 
a high binary frequency: near-infrared adaptive optics studies by 
Chauvin et al. (2002 ) and Brandeker et al. (2003 ) revealed that 
MBM 12 contains at least 4 binaries and 2 triples with projected 
separations between 20 and 4000 AU (of which one is a can- 
didate quadruple: LkHa 262 has a projected distance of only 
15 arcsec to the triple LkHa 263). 

In this paper, we present Spitzer IRS and MIPS observations 
of all MBM 12 members, and combine our observations with lit- 
erature photometry and spectroscopy to analyse the derived dust 
characteristics by relating them to the stellar parameters and disc 
properties. In Sect. [2] we present the individual targets and the 
Spitzer observations. The analysis in Sect. [3] first discusses the 
spectral energy distributions, and then the spectral features are 
interpreted in terms of dust properties with the aid of a spectral 
decomposition model. In Sect. [4] we discuss grain growth, crys- 
tallisation, disc properties and accretion rates. We round off with 
conclusions in Sect. [5] and give more details in the Appendix on 
the dust model used and the mass fractions derived. 

2. Sample and Observations 

2.1. Targets 

Our sample is unbiased, as it includes all the confirmed members 
of the MBM 12 cloud (Luhman I200T1 [2004b . Their parameters, 
spectral type and T Tauri Class (weak-line or classical) are listed 
in Table Q] We re-derived their T Tauri Class, following White 
and Basri (2003): T Tauri stars are classical when, for spectral 
types between K0 and K5, the equivalent width of their Ha line, 
|EW(Ha)| > 3 A, for K7 to M2.5, |EW(Ha)| > 10 A and for 
M3 to M5.5, |EW(Ha)| > 20 A. We also list their multiplicity 
status and projected separation of companions (when present), as 
determined by Brandeker et al. (2003 ) and Chauvin et al. (2002 ) 
for 8 of the 12 members. The objects under consideration are 
all T Tauri stars with spectral types between K3 and M5.75, and 
at least 6 out of the 12 targets are known to have companions. 
Below, we list more information on the individual targets: 

MBM 12-1 (RX J0255.4+2005) is a weak-line binary, with 
the primary of spectral type K6 and a separation of 0."533 
(Chauvin et al. [2002b . X-ray observations with ROSAT 



showed that this object flares, with a factor six increase in X- 
ray counts during the flare. Spectral fits of the count rate sug- 
gest a coronal origin, not untypical for other flaring WTTS 
(Hearty et al. |2"000j ). 

LkHa 262 is a CTTS, with a disc detected at 2.1 mm, from 
which Itoh et al. (120031 ) estimated a disc mass of 0.05 M . It 
is located at a distance of only 15. "3 from LkHa 263. 

LkHo- 263 is a triple system, with LkHa 262 possible belong- 
ing to this triple to then form a quadruple system (Chauvin 
et al. 2002); however, it is not clear whether this system 
is bound. The C component of LkHa 263 has spectral 
type M0, and was found to harbour an optically thick disc, 
that spatially-resolved observations showed to be edge-on 
(Jayawardhana et al. 120021 ). These authors derived a disc 
mass of 0.0018 M for this 0.7 M C component, based on 
their near-IR images. Furthermore, forbidden lines in the op- 
tical spectrum of LkHa 263 C suggest the presence of a jet 
(Jayawardhana et al. l2002b . 

LkHa 264 is a wide binary with a separation of 9." 160. 
Millimetre observations at 1.3 and 2.1 mm give a disc mass 
of 0.09 M around the primary (Itoh et al. 2003). Emission 
of molecular hydrogen at 2. 1218 and 2.2233 pm was also de- 
tected around the primary. The width of these lines points to 
a disc origin, while further modelling locates the NIR emit- 
ting H2 in the inner 10 AU, and shows that the disc is seen 
nearly pole-on (Carmona et al. 120081) . The line strength ra- 
tio is consistent with a temperature lower than 1500 K, and 
points to thermal excitation by UV photons; LkHa 264 has 
a strong UV excess, so it is indeed plausible that there are 
enough UV photons to excite the H2. The total mass of the 
optically thin, hot H2 in the disc of the primary is estimated 
to be a few lunar masses (Carmona et al. 2008 ). 

MBM 12-5 (E 0255+2018) is a binary TTS with a separation 
of l."144. Following the classification by White & Basri 
(2003 ), its spectral type of K3 and Ha equivalent width of 
3.1 A puts it on the border between the classical and weak- 
line T Tauri stars. 

MBM 12-6 is a classical T Tauri star for which Brandeker et 
al. ([2003 ) did not find a companion within a radius of l."6. 

MBM 12-7, 9 and 11 are all weak-line T Tauri stars, with 
the latest spectral types of the whole MBM 12 sample: be- 
tween M5 and M6, corresponding to masses between 0.15 
and 0. 1 M Q . We could not find any multiplicity data for these 
relatively faint objects. High spatial resolution images are 
needed to establish the multiplicity status of these sources. 

MBM 12-8 is a classical TTS, with spectral type M5.5, and 
the highest Ha equivalent width: 120 A, suggesting it is the 
most actively accreting object of the sample. No additional 
information concerning its multiplicity status is available. 

MBM 12-10 is a binary weak-line TTS with a separation of 
0."390. 

S 18 is a triple system, consisting of the primary A and, at 
a projected distance of 0."747, a tight binary companion 
Bab (with a separation of only 0."063 - Brandeker et al. 
2003 ). Millimetre observations at 2.1 mm give a disc mass 
of 0.07 M G around the primary (Itoh et al. |2003b . 



2.2. Spitzer IRS spectroscopy 

The T Tauri stars in MBM 12 were observed with Spitzer as 
part of a larger programme on young stellar clusters to study 
the evolution of protoplanetary discs (GO proposal 20691, RI. 
Bouwman). We obtained 7.5 - 35 pm low-resolution (R = 
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Table 1. Target coordinates and parameters: Spectral Type, T e g, Ha equivalent width of the primary from Luhman (2001), T Tauri 
Class derived by us (see Sect. 12. IK multiplicity (S: single, B: binary, T: triple and Q: quadruple) and projected separations from 
Chauvin et al. (2002 ) and Brandeker et al. (2003). A negative equivalent width means emission; the width is given in angstroms, A. 



Object 


Alternative Name 


a-(2000.0) 


(5(2000.0) 


Spectral 


T c ff 


Eq. Width 


Class 


Multiplicity Projected 






(h m s) 


C ' ") 


Type 


(K) 


(Ha) 




Status Separation (") 



MBM 12-1 


RXJ0255. 4+2005 


02 55 25.78 


20 04 51.7 


K6 


4205 


-1 


WTTS 


B 


0.533 


MBM 12-2 


LkHor 262 


02 56 07.99 


20 03 24.3 


M0 


3850 


-40 


CTTS 


(OX 


15.3 


MBM 12-3 


LkHa 263 (A-B/B-C) 


02 56 08.42 


20 03 38.6 


M3 


3415 


-25 


CTTS 


T 


0.416/4.1 


MBM 12-4 


LkHa 264 


02 56 37.56 


20 05 37.1 


K5 


4350 


-18 


CTTS 


B 


9.160 


MBM 12-5 


E 02553+2018 


02 58 11.23 


20 30 03.5 


K3 


4660 


-3 


W/CTTS* 


B 


1.144 


MBM 12-6 


RX J0258.3+1947 


02 58 16.09 


19 47 19.6 


M5 


3200 


-29 


CTTS~ 


S 




MBM 12-7 


RX J0256.3+2005 


02 56 17.98 


20 06 09.9 


M5.75 


3024 


-14 


WTTS 






MBM 12-8 




02 57 49.02 


20 36 07.8 


M5.5 


3058 


-120 


CTTS 






MBM 12-9 




02 58 13.37 


20 08 25.0 


M5.75 


3024 


-10 


WTTS 






MBM 12-10 




02 58 21.10 


20 32 52.7 


M3.25 


3379 


-12 


WTTS 


B 


0.390 


MBM 12-11 




02 58 43.80 


19 40 38.3 


M5.5 


3058 


-14 


WTTS 






MBM 12-12 


S 18 (A-B/Ba-Bb) 


03 02 21.05 


17 10 34.2 


M3 


3415 


-69 


CTTS 


T 


0.747/0.063 



" Possible companion to LkHa 263 ABC 

* On the border between weak-line and classical TTS, see target notes in Sect. 12. II 



60 - 120) spectra of the MBM 12 cluster members with the 
Infrared Spectrograph (IRS, Houck et al. 2004) on-board the 
Spitzer Space Telescope. A high accuracy PCRS peak-up was 
executed prior to the spectroscopic observations to position the 
target within the slit. All targets have been observed with a 
minimum of three observing cycles for redundancy. Our spec- 
tra are based on the droopres products processed through the 
S 15.3.0 version of the Spitzer data pipeline. Partially based on 
the SMART software package (Higdon et al. 2004), our data was 
further processed using spectral extraction tools developed for 
the "Formation and Evolution of Planetary Systems" (FEPS) 
Spitzer science legacy team (see also Bouwman et al. [2008). 
The spectra were extracted using a 6.0 pixel and 5.0 pixel fixed- 
width aperture in the spatial dimension for the observations with 
the first order of the short- (7.5-14 pm) and the long- wavelength 
(14-35 pm) modules, respectively. The background was sub- 
tracted using associated pairs of imaged spectra from the two 
nodded positions along the slit, also eliminating stray light con- 
tamination and anomalous dark currents. Pixels flagged by the 
data pipeline as being "bad" were replaced with a value inter- 
polated from an 8 pixel perimeter surrounding the errant pixel. 
The low-level fringing at wavelengths > 20 pm was removed 
using the irsfringe package (Lahuis & Boogert 2003). To re- 
move any effect of pointing offsets perpendicular to the slit, we 
matched orders based on the point spread function of the IRS 
instrument, correcting for possible flux losses (see Swain et al. 
2008 for further details). To remove any effect of pointing offsets 
parallel to the slit, we made a flatfield correction to the nominal 
rsrf which was derived from calibration measurements of stan- 
dard stars where the calibrator was mapped along the slit. The 
spectra are calibrated using a spectral response function derived 
from multiple IRS spectra of the calibration star rji Doradus and 
a MARCS stellar model provided by the Spitzer Science Centre. 
The spectra of the calibration target have been extracted in the 
same way as our science targets. The relative errors between 
spectral points within one order are dominated by the noise on 
each individual point and not by the calibration. We estimate a 
relative flux calibration across an order of « 1 % and an abso- 
lute calibration error between orders/modules of « 3 %, which is 
mainly due to uncertainties in the scaling of the MARCS model. 



Table 2. MIPS photometry at 24 and 70 pm, and their statistical 
errors (lcr flux uncertainties, derived as described in the text). 
For those objects that were not detected, 3cr upper limits were 
derived (see Sect. 12. 31 ). The absolute calibration uncertainties are 
4 and 7% for 24 and 70 pm, respectively. 



Object 


24 ytim (lcr) 


70 yum(lcr) 




(mJy) 


(mJy) 


MBM 12-1 


2.8 (0.4) 


< 34.4 


MBM 12-2 


142 (0.4) 


216° (12) 


MBM 12-3 


50.8 (0.4) 


216~(12) 


MBM 12-4 


282 (0.5) 


266 (13) 


MBM 12-5 


308 (0.5) 


253 (12) 


MBM 12-6 


25.7 (0.4) 


< 34.4 


MBM 12-7 


< 1.2 


<32.9 


MBM 12-8 


< 1.2 


<41.4 


MBM 12-9 


< 1.2 


<37.5 


MBM 12-10 


21 (0.4) 


<35.6 


MBM 12-11 


< 1.2 


< 34.1 


MBM 12-12 


45.5 (0.4) 


< 32.2 



" Combined flux of LkHo- 262 and 263 



The targets with the latest spectral types, MBM 12-7, 8, 9 
and 11, were not detected with IRS, so we will not include them 
in the rest of our analysis of the dust properties. This means that 
the spectral range for which we have IR spectroscopy is limited 
to K3-M5. 

2.3. Spitzer MIPS photometry 

Additional infrared data was obtained using MIPS (Rieke et al. 
2004) on Spitzer (Werner et al. 120041) in photometry mode at 24 
and 70 pm (default scale). All of our targets were observed in 
February 2006, using 1 cycle of 3 sec integrations at 24 pm and 
2-5 cycles of 10 sec integrations at 70 pm, corresponding to on- 
source intergation times of 24. 1 sec and 25 1 .6 - 629 sec at 24 and 
70 pm, respectively. Our observations were processed using ver- 
sion S 16.0.1 of the Spitzer Science Center (SSC) data pipeline. 
We created two 24 pm mosaics for each of our sources from the 
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resulting basic calibrated data (BCD) images of each data collec- 
tion event ( DCE) using the SSC's MOPEX software (Makovitz 
& Marleau 120051 1. one with a pixel scale approximately that of 
the native pixel scale (2. "45 pixel -1 ) and another resampled to 
approximately one half of the native pixel scale (1 ."23 pixel -1 ). 
We, similarly, created two 70 pm mosaics for each of our sources 
using the filtered BCDs (which have spatial and temporal filters 
applied to the data in order to remove instrumental signatures) 
with outlier rejection, one with a pixel approximately that of the 
native pixel scale (9. "9 pixel -1 ) and another resampled to ap- 
proximately one quarter of the native pixel scale (2."0 pixel -1 ). 

We used the APEX portion of MOPEX to perform aperture 
photometry on our mosaicked images. APEX applies a median 
filter to the data to estimate the sky background at any pixel in 
the image and subtracts the median-filtered image before sum- 
ming the flux in an aperture. Since the estimated backgrounds 
in the majority of our fields, extrapolated from COBE (Cosmic 
Microwave Background Explorer), was medium to high, we used 
a medium-sized aperture with radius of 6" at 24 pm and 16" at 
70 /im to reduce the amount of background contamination in the 
aperture. These apertures are not large enough to contain all of 
the photons from a diffraction-limited point source; therefore, 
we applied scalar aperture corrections of 1.697 and 1.771 at 24 
pm and 70 /im, as published on the SSC website. The 70 /im 
aperture correction is somewhat dependent on the colour of the 
source; the aperture correction used here assumes that the mea- 
sured flux has a red power law shape, F v cc y -2 , because any 
flux detected at 70 /im is dominated by the emission from a 
cool, dusty disc (T ~ 100 K). The MIPS data handbook (ver- 
sion 3.2.1) states that products processed with version S14.4 of 
the data pipeline have absolute calibration uncertainties of 4 and 
7% at 24 and 70 //m, respectively. 

The majority of our sources were detected with signal to 
noise ratios greater than 10 at 24 pm and were not detected at 
70 pm (see Table [2J. We estimate lcr flux uncertainties for ob- 
jects that were detected and 3cr upper limits on the fluxes of ob- 
jects that were not detected from uncertainty mosaics produced 
by MOPEX (at the native pixel resolution). To determine the lcr 
uncertainty in a measured flux, we take the square root of the 
sum of the uncertainties (in the flux in each pixel in the aperture) 
in quadrature, multiplied by the aperture correction (1.164 and 
1.197 at 24 or 70 pm, respectively), centered at the expected po- 
sition of the source. Similarly, we determine the 3cr upper limit 
on the 24 and 70 pm fluxes as three times the square root of the 
sum of uncertainties (in the flux in each pixel in the aperture) 
in quadrature, multiplied by the aperture correction (1.164 and 
1.197, respectively), in a relatively large aperture, with radius 
35", centered at the expected position of the source. 

In Table |H we list the MIPS photometry for the MBM 12 
members. The 4 stars with the latest spectra types, MBM 12-7, 
8, 9 and 1 1 were not detected at 24 nor at 70 //m. In Fig. [2] we 
show that the MIPS photometry at 24 /im and the fluxes from the 
IRS spectra agree very well. 

3. Analysis 

3. 1 . Spectral energy distributions 

To complement our Spitzer observations, we collected photom- 
etry from the literature: Jayawardhana et al. (2001 ) and Luhman 
d2U0Tb for RI, JU K S and LN band photometry. The V band is 
from Broeg et al. (2006), while the submillimetre and millime- 
tre photometry are from Hogerheijde et al. (120031 1 and Itoh et 
al. ([2003 ), respectively. In Fig. [2] we show the spectral energy 



2.5 
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0.0 i i i 

10 20 30 40 

X [/im] 

Fig. 1. IRS Spitzer spectra of the seven MBM 12 sources that 
show excess emission in the infrared. The strength and shape of 
the emission features, as well the slope of the spectra varies from 
source to source. 



distributions (SED) of the 12 targets, together with a MARCS 
stellar model (Gustafsson et al. 2008 ) and for those sources with 
r e ff> 4200 K with a Kurucz atmosphere model (Kurucz 1994) 
for the appropriate effective temperature of the central source, as 
listed in Table Q] The photometry was dereddened using the Ay 
values derived by Luhman (2001), assuming a standard extinc- 
tion law (Ry = 3.1). 

Of the 8 sources detected with Spitzer, only one object, 
MBM 12-1, a K6-type weak-line TTS with a companion at 
0."533, shows no excess emission at all: the spectrum is purely 
atmospheric. The other 7 sources show an infrared excess, and 
also have a silicate emission feature at 10 pm, with varying de- 
grees of strengths and shapes (see Fig. [TJ. As we do not have 
any data between ~3.4 and 8 pm, it is difficult to determine the 
properties of the inner disc, e.g. whether there is an inner hole 
or not, as witnessed by a lack of excess shortwards of the 10 mi- 
cron feature. Sources MBM 12-3, 6 and 10, appear not to have a 
near-IR excess (see Fig . [2j ; they are candidate transitional discs, 
where the inner disc has already been cleared out. MBM 12-2, 
4, 5 and 12 have 'full discs', as they already have an excess in 
the K band. Additional data, between 3 and 8 micron, is needed 
to confirm the status of the transitional candidates. An overview 
of the main features observed is given in Table [3j the properties 
of the dust, as derived from the solid-state features, will be dis- 
cussed in the following section. We discuss the relation between 
excess luminosity and binarity in Sect. 14.31 

3.2. Dust mineralogy 

The most frequently analysed region in the context of dust prop- 
erties is the 10 micron region: it is through this window that dust 
characteristics like composition (e.g. olivines, Mg2 X Fe2(i- X )Si04 
and pyroxenes, Mg x Fe(i_ x )Si03, or Mg/Fe content), structure 
(crystalline vs. amorphous) and grain sizes can be derived from 
their solid-state features. The 10 pm region is only sensitive to 
dust grains up to sizes of a few microns, as the feature flattens 
out when the dust grains attain sizes similar to the wavelength, 
making them essentially invisible. Furthermore, it is important 
to realise that the dust causing this feature traces only a small 
fraction of the disc material, namely those dust grains that are 
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Fig. 2. The spectral energy distributions of the 12 T Tauri stars in MBM 12. The lowest-mass stars, MBM 12-7, 8, 9 and 11 were 
not detected with Spitzer, therefore we can only show upper limits in the infrared for those sources. MBM 12-1 is the only source 
in our detected sample that shows no excess above its atmosphere. MBM 12-3, 6 and 10 appear to not have a near-IR excess, so are 
candidate transitional discs, while MBM 12-2, 4, 5 and 12 have 'normal accretion discs', with an excess extending into the near-IR. 



located in the optically thin disc atmosphere, while the bulk of 
the dust mass is located in the disc midplane. In addition, as the 
dust temperature needs to be high enough (150-450 K) in order 
to radiate at 10/im, the radial location of the dust observed is also 
limited: for a typical star in our sample, this is around 1 AU. At 
longer wavelengths, features of crystalline dust (forsterite and 
enstatite) witness another important part of dust processing in 
the protoplanetary disc. 

In Fig. [3] we zoom in on the three most interesting ranges 
that are observed with IRS, and also indicate the location of the 
most important features. Apart from the broad feature of amor- 
phous silicate at 9.7 (present in all 7 stars with a detected ex- 
cess) and 18 /mi, features of enstatite (MgSiOs) and forsterite 
(Mg2SiC>4) are observed. We do not see evidence for carbona- 
ceous dust, such as features from polycyclic aromatic hydrocar- 
bons (PAHs). This is perhaps not surprising, given that PAHs are 
transiently excited by UV photons and the central sources have 
low temperatures, although some sources do show a rather large 
UV excess (e.g. MBM 12-2, 4 and 12). 



In order to analyse the composition of the dust in the disc 
atmosphere, the radiation of which dominates the IRS spectrum, 
we use the two-layer temperature distribution (TLTD) spectral 
decomposition routines described in Juhasz et al. ((2009). This 
method uses a multi-component continuum (star, inner rim, disc 
midplane), assuming that the region where the observed radia- 
tion originates (both optically thin and thick) has a distribution 
of temperatures. In this fitting method, the observed flux-density 
at a given frequency is given by 

N M r T "- M 2n 

F v = F v , cont + YjYjDijKijj ' -^B v (T)T^dT (1) 

where N and M are the number of dust species and of grain 
sizes used, respectively, kij is the mass absorption coefficient of 
the dust species i and grain size j. B V (T) is the Planck-function, 
qa is the power exponent of the temperature distribution and d 
is the distance to the source. The subscript 'a in the integration 
boundaries refers to the disc atmosphere. The continuum emis- 
sion (F V) cont) is given by 
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Fig. 3. Zoom in on the 3 main ranges in the IRS spectra: 7-14 /mi, 17-25 /mi and 26-36 /mi. We overplot the positions of the main 
crystalline features: long-dashed lines indicate enstatite, while dot-dashed lines indicate forsterite. 



^v.cont = Dq^BAT*) + Dl 



+ D2 



r TrM 2n 
J, 



B v (T)T~dT (2) 
B v (T)T^dT. (3) 



The first term on the right hand side describes the emission 
of the star, while the second and third term describe the radia- 
tion of the inner rim (subscript V) and the disc midplane (sub- 
script ';«'), respectively. The meaning of the different parameters 
are summarised in Table [7] in the Appendix. For each compo- 



nent (disc atmosphere, inner rim, midplane), the highest temper- 
ature is fitted to obtain T a / r / m> - m , while the lowest temperature, 
Ta/r/m, out, is calculated requiring that the annulus with that tem- 
perature contributes more than 0. 1 % to the total flux. 

For our fit, we used five dust species that are commonly 
found in the discs of young stars, in three grain sizes (0.1, 1.5 and 
6 /mi). We applied the theory of distribution of hollow spheres 
for the crystalline dust (Min et al. 120051 and the classical Mie 
theory for spherical particles for the amorphous dust, to derive 
the mass absorption coefficients from the optical constants. The 
list of the dust species, the origin of the optical constants and the 
grain model used are presented in Table [4] Furthermore, since 
the dust grains radiating at 8-13 /mi and 20-30 /mi are at dif- 
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Fig. 4. The TLTD fits of the 10 micron silicate feature (thick, solid line), and the IRS spectra with their noise (black points with error 
bars). Notice the different shapes that are present in this small sample: MBM 12-5 has a triangular shape, pointing to small, amor- 
phous silicate grains, while MBM 12-6 has a much broader feature, indicating larger and crystalline grains. The narrow emission 
feature at 9.3 fim is caused by crystalline enstatite. 



Table 3. Summary of the spectral appearance of the detected tar- 
gets. We list the shape of the spectral energy distribution, where 
we distinguish between a 'full' and a 'transitional' disc, based 
on the presence or absence of excess emission at wavelengths 
shorter than 8 micron, and write 'transitional object (TO) can- 
didate' when this classification needs to be confirmed by addi- 
tional photometry. Furthermore, we list the wavelengths at which 
features of enstatite and forsterite are present. Colons indicate 
uncertain detections. 



Object 


SED 


Enstatite 


Forsterite 


MBM 12-1 


Photosphere 






MBM 12-2 


Full 




19, 27.8, 33 fim 


MBM 12-3 


TO candidate 


9.3: 


11.3:, 19:, 27.8, 33 fim 


MBM 12-4 


Full 




11.3:, 19:, 27.8, 33 fim 


MBM 12-5 


Full 




19, 27.8, 33 fim 


MBM 12-6 


TO candidate 


9.3 yum 


27.8, 33 fim 


MBM 12-10 


TO candidate 


9.3 fim 


23.4, 27.8 fim 


MBM 12-12 


Full 


9.3 jum 


11.3, 19, 27.8, 33 fim 



Table 4. Overview of dust species used in our fitting routines. 
For each component, we specify its lattice structure, chemical 
composition, shape and reference to the laboratory measure- 
ments of the optical constants. For the homogeneous spheres, 
we used Mie theory to calculate the opacities. For the inhomoge- 
neous spheres, we used the distribution of hollow spheres (Min 
et al. 120051) . to simulate grain shapes deviating from perfect sym- 
metry. 



Species 



State Chemical 



Shape 



Amorphous silicate 1 
(Olivine stoichiometry) 
Amorphous silicate' 
(Pyroxene stoichiometry) 
Forsterite 2 
Clino Enstatite 3 
Silica 4 



formula 


A 


MgFeSi0 4 


Homogeneous 






Sphere 


A 


MgFeSi 2 6 


Homogeneous 






Sphere 


C 


Mg 2 Si0 4 


Hollow Sphere 


C 


MgSi0 3 


Hollow Sphere 


A 


Si0 2 


Hollow Sphere 



References: (1 ) Dors chner et al. 119951 (2) Se rvoin et al. [19731 
(3) Jager et al. [19981 (4) Henning et al. fTWTI 



ferent temperatures, hence at different radial distances from the 
star, we splitted the wavelength region in two parts before fitting. 
This way, we can take into account that different species can 
contribute in different amounts to the spectrum, if there is a ra- 
dial gradient in their abundances and/or properties. We chose the 
1) the 7-17 fim range for the shorter wavelength region, which 
Juhasz et al. ((2009) showed to be the optimal wavelength range 
to balance the presence of spectral features and radial changes, 
and 2) the 17-37 fim range for the longer wavelength region, 
which includes important crystalline features; both regions were 
separatedly fitted, and we will refer to them as the 'shorter' and 
'longer' wavelength regions. 



As the IRS spectra do contain noise, we repeated the fit of 
each object 100 times, every time adding random Gaussian noise 
on top of the measured spectrum. The final composition is de- 
rived as the average of those 100 fits, while the errors are derived 
from the standard deviation, taking into account the positive and 
negative directions of the variations. Therefore, the errors we 
give describe the S/N and measurement errors in the spectra. We 
emphasize that we cannot give error estimates that are related to 
species that are not included in our fits, nor of different shapes 



8 



G. Meeus et al.: MBM 12: young protoplanetary discs at high galactic latitude 



Table 5. Distilled fitting results: given are mass-averaged size 
of the amorphous and the crystalline silicates, < a > a m.sii. and 
< a > cr yst.sii., respectively, mass fraction of the crystalline grains, 
fcryst (ratio of crystalline silicate to total silicate mass), and the 
mass ratio of forsterite to enstatite, as derived for both the 7-17 
and 17-37 pm regions. 



Object 


< a > am .sil. 


< a >cryst.sil. 


f 

Hrvsl. 


forst/enst 




(jim) 


Oum) 


(%) 





MBM 12-2 
MBM 12-3 
MBM 12-4 
MBM 12-5 
MBM 12-6 
MBM 12-10 
MBM 12-12 



4-1 O 

4-0 dp 

0.4 O 



-17 pm region 
5-2 O 

1.7 dl) 



2.2 C 
3.9 C\ 



5.3 (!) 
6.0 (!« *) 



4-1 C : ) 
4.1 (!) 
5.3 Cf ? ) 
1.9 CS) 
5.5 (i{) 



10-7 Q 

16.3 C^) 
2.9 

46.5 Cp 
4.6 O 
32.8 (!") 



o.i O 

0.4 H») 

o.i dto 

0.3 Cjj) 

o.o cj) 

0.8 (!) 

o-i ds) 







17-37 pm region 






MBM 12-2 


o-i O 


o.8 CSS) 


1-8 O 




MBM 12-3 


1.2 d) 
0.2 d \) 


o.4 db 


4.8 dl) 


1.2 dp 


MBM 12-4 


1-6 C |) 


2.8 dp 


i.o d) 


MBM 12-5 


o.7 d?) 


4.0 dj) 


3.5 (!>S) 


0.9 


MBM 12-6 


i.5 dp 


o-i C|) 


7.1 c;1) 




MBM 12-10 


0.4 (! ) 




2.9 cfi) 
3.3 (!>*) 




MBM 12-12 


0.6 (! ) 


0.8 CJ«) 


0.9 O 



or composition of the grains. Those errors are subject of another 
study (A. Juhasz et al. 2009, in preparation). 

With the fits, we determine the source of the solid-state fea- 
tures that are observed in the spectra. Since several species (e.g. 
carbon) are featureless within our wavelength range (but can 
contribute to the continuum), they remain undetected in our fits; 
so it is important to realise that we can only discuss the 'visible 
dust', that has features in the range 7 to 37 micron. In Fig [4] we 
show the best fit of the 7 to 17 micron region, and a summary of 
the derived dust properties for both the shorter and longer wave- 
length regions is given in Table [5] The detailed results of the fits 
in both regions, in terms of mass fractions of the relevant species 
and their sizes, are given in Table[8]in the Appendix. 

4. Discussion: 2 Myr-old T Tauri discs 

4. 1 . Grain growth 

A first thing to notice, when looking at the 10 micron silicate 
feature, is that the strength of the emission varies. Van Boekel et 
al. ((2003) showed that the shape and the strength of the 10 pm 
feature in Herbig Ae/Be stars are related, and demonstrated it 
to be evidence for grain growth: a strong and triangular 10 pm 
feature is typical of small (submicron-sized) grains, whereas a 
weaker and broader feature indicates the presence of larger sized 
grains (up to a few microns). In Fig. [5] we probed this relation 
for our sample of T Tauri stars, and also found evidence for grain 
growth in the protoplanetary discs of the MBM 12 members. 

Kessler-Silacci et al. ( 120061120071) and Apai et al. d2005l) no- 
ticed that later type TTS have weaker 10 micron features than 
earlier type TTS, and related this to the location of the dust caus- 
ing this feature. This relation was confirmed by Sicilia-Aguilar 
et al. ( 120071 ), who found that the presence of a very weak feature 
was 3 times more frequent for M-type stars than for earlier-type 
stars. The temperature of the silicates causing the feature is ~ 
300 K, a temperature that is reached more inwards for lower lu- 



I .2 - 



0.8 



0.6 - 



1.0 1.5 2.0 2.5 3.0 

Peak over continuum (10 micron feature) 



3.5 



Fig. 5. Relation between the shape of the 10 pm feature (flux 
ratio 11.3 over 9.8 pm) and the feature strength (peak over con- 
tinuum). Our objects (squared symbols, numbered as listed in 
Table[T]l, follow the trend observed for TTS in Tr 37 (4 Myr) and 
NGC 7160 (12 Myr), shown by diamonds (Sicilia-Aguilar et al. 
120071 ) and TTS in Taurus (2 Myr), shown by triangles (Pascucci 
et al. 120081 ): the flux ratio is correlated with the feature strength. 



minosity sources than for higher luminosity sources: under the 
assumption of a similar disc structure for all spectral types dis- 
cussed here, an MO-type star reaches 300 K between 0.7 and 
1.5 AU, while it is between 1.2 and 3 AU for a K5-type star, a 
factor 2 difference in distance. Furthermore, the density distri- 
bution has a radial dependence, as it decreases with increasing 
radius, and grain growth occurs faster in more dense environ- 
ments. This implies that grain growth will be faster in the more 
inward regions, so that lower luminosity sources will appear to 
have larger grains, as their 10 micron feature originates from 
a denser region, in which growth naturally occurs more rapidly. 
This is also seen in our spectra: in Fig. [6] we plot the continuum- 
normalised spectra, and also list their spectral types. With the 
exception of MBM 12-4, the fainter objects indeed have weaker 
features. This is further quantified through our fitting with the 
TLTD method: in Fig. [7] we show the relation between the mass- 
averaged grains size in the 10 micron region and the effective 
temperature of the central star. There is a clear trend for cooler 
stars to have larger grain sizes. 

Of course, not all discs need to have the same structure, and 
the degree of flaring can also play an important role in this con- 
text. However, to study these effects, detailed radiative transfer 
modelling of each source is needed, which is beyond the scope 
of this paper. 

4.2. Crystallisation 

The dust that is initially incorporated into the protoplanetary disc 
is largely amorphous, as it comes from the ISM, for which an 
upper limit of ~2% in mass of the crystalline grains was deter- 
mined using spherical grains (Kemper et al. 120041) . Min et al. 
(120071 ) further improved this number using grains with irregu- 
lar shapes and derived a crystallinity of only However, the 
amorphous dust may become crystalline through thermal anneal- 
ing (e.g. Fabian et al. 120001 ) or shock heating (e.g. Scott & Krot 
2005 ). We derived the crystalline mass fraction in the shorter 
wavelength region, and found a large degree of variation: be- 
tween 2.9!| 2 % for MBM 12-5 and 47+^ % for MBM 12-6. 
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Fig. 8. MBM 12-6, host of a large amount of crystalline silicates, 
witnessed by the agreement in peak positions between the spec- 
trum and the features of the mass absorption coefficients of en- 
statite (Jager et al. 1998 ; dotted line) and forsterite (Servoin et al. 
119731 dashed line). For comparison, we also show the Herbig Be 
star HD100546 (Malfait et al. |1998l l, and the solar-system comet 
Hale-Bopp (Crovisier et al. 1997 ), both host of large amounts 
of crystalline dust. The fluxes of HD 100546 and Hale-Bopp are 
scaled to match the flux of MBM 12-6. 
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Fig. 6. Continuum normalised flux in the 10 fim region. The ob- 
jects with the latest spectral types tend to have the weakest fea- 
ture. 
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Fig. 7. Mass-averaged size of the amorphous silicate grains in 
the warmer region, versus effective temperature of the central 
star. The cooler the star (the later the spectral type), the larger 
the dust grains tend to be. 



The crystalline fraction of MBM 12-6 is not very well deter- 
mined, therefore we use the lower limit of 30 % that was found 
for this object. The mass fraction of the crystalline silicates is 
not related to the spectral type. Furthermore, we do not find a 
correlation between crystallinity and the size of the amorphous 
grains; also the size of the crystalline grains appears to be un- 



correlated with amorphous grain size (see Fig. [9j>. However, we 
do note that the sources with the largest crystalline mass frac- 
tions (MBM 12-6 and MBM 12-12) also have the largest grains 
size, both for the amorphous as for the crystalline grains (see also 
Fig. [Toll. In the longer wavelength region, we find the variation in 
crystallinity to be less, and the crystalline mass fraction smaller: 
between 1.8+JJ | % for MBM 12-2 and 7.1+^ % for MBM 12- 
6. This difference in crystallinity between both regions for the 
whole sample could be related to the presence of a large amount 
of small amorphous silicate grains in the cooler disc region (av- 
erage size for the sample is 0.7 /mi), while in the warmer disc 
region the amorphous silicate grains are much larger (average 
size for the sample is 3.7 //m). Larger grains will produce less 
detectable features, so that the crystalline grains that still show 
features will appear more abundant. 

In Fig. [8] we show the 15 to 37 fim spectrum of MBM 12-6, 
the object with the largest fraction of crystalline grains. Its IR 
spectral appearance is similar to that of HD 100546, a Herbig 
B9e star with highly evolved dust, which is remarkable, given 
that their effective temperatures are very different: 1 1000 K for 
the B9 star, and only 3200 K for the M5 type star. Still, the dust 
around both objects is dominated by crystalline grains, indicat- 
ing that the temperature of the central star does not play an im- 
portant role in the crystallisation process. The spectrum of MBM 
12-6 is also similar to that of solar-system comet Hale-Bopp, so 
that similar dust processing mechanisms must exist in both our 
solar system and the disc of MBM 12-6. 

In the dust model by Gail (1998), the composition of the dust 
is derived, based on condensation sequences and chemical equi- 
librium considerations. He predicts that - assuming crystalline 
silicates form as high temperature gas phase condensates - of the 
crystalline silicates, enstatite will be the dominant constituent, 
while forsterite is only present in a small region close to the star. 
More recent simulations show that radial mixing within the disc 
is capable of transporting a substantial amount of crystalline ma- 
terial from the inner disc towards more colder, outward regions 
(Keller & Gail 120041 
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Fig. 9. Relation between the properties of the crystalline silicates 
and the size of the amorphous silicates in the 7 to 17 micron re- 
gion. Upper panel: mass-averaged size of the crystalline silicates 
versus mass-averaged size of the amorphous silicates. The two 
appear unrelated. Lower panel: crystalline mass fraction versus 
mass-averaged size of the amorphous silicates. Also here we do 
not see a correlation. 
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Fig. 10. Properties of the the crystalline silicates, derived from 
the 7 to 17 micron region. Upper panel: the forsterite over en- 
statite ratio in function of the size of the crystalline silicates. 
There is a trend for larger grains to have a smaller forst/enst ra- 
tio. Lower panel: crystallinity in function of crystalline grain size 
. The crystalline fraction only exceeds 20% when the crystalline 
grains are larger than 5 micron. 



To better understand the crystallisation process occurring in 
the T Tauri discs, we look at the forsterite to enstatite ratio. In 
Table [5] we list the forsterite to enstatite mass ratio for both re- 
gions. In the warmer region, the smallest ratio is found for those 
sources with the largest grains (see Fig. [TOt . implying that en- 
statite is more abundant (relative to forsterite) when larger grains 
are present. We further compare the relative mass fractions of 
forsterite and enstatite in both regions: the average ratio for the 
sample is 0.3 in the warmer region, and 0.9 in the cooler region. 

It is thus clear that there is a spatial gradient in the forsterite 
to enstatite ratio, with forsterite dominating the cooler regions, 
and enstatite more abundant in the warmer, inner regions. Such 
a gradient was already found in a few earlier studies, e.g. 
Bouwman et al. ( 120081 ). This is in contrast with the predictions 
by Gail (1998 2004), suggesting that the chemical equilibrium 
conditions needed for the forsterite to enstatite conversion are 
not reached inside these discs, and that the crystallisation process 
must be different, even when including radial mixing. It is inter- 
esting to consider that enstatite is more abundant in the warmer 
(inner) regions, and that we also find a larger enstatite abundance 
in those sources that have larger grains (see Fig. ITOt. These ob- 
servations suggest that enstatite forms easier in regions of higher 
density, where also grain growth is more abundant. Bouwman et 
al. ([2008 > discuss the formation of enstatite and forsterite in de- 
tail, and link the formation of enstatite in the inner region with 
the conditions that prevail there: due to the higher density and 
temperature, it takes the dust grains longer to cool down so that, 
potentially, equilibrium conditions can be reached - what is more 



unlikely at larger radial distances, where both the density and the 
temperature are lower. 

4.3. Disc properties 
4.3.1 . Disc fraction 

We obtained Spitzer data for the complete sample of MBM 12 
members; 4 of them, however, could not be detected. The re- 
maining 8 objects (spectral type range K3 to M5) have a disc 
fraction rate of 7 out of 8, or nearly 90%. This is very high, 
when compared to a similar spectral type range in other star 
forming regions: Damjanov et al. ( 120071 ) derive a disc fraction 
of 47% for M0-M4 and 55% for K3-K8 type stars in the 2 
Myr-old Chamaeleon I and 63% region; however, Flaherty and 
Muzzerole (2008 ) derive a disc fraction between 65 and 81% for 
the 2 Myr old clusters NGC 2068 and NGC2071. The high disc 
fraction rate found for MBM 12 may (partly) be attributed to the 
absence of close companions to the disc bearing stars: there are 
no companions found at a projected distance smaller than 0."39 
(or 70 AU at the distance of MBM 12). 

Of the 7 T Tauri stars with a disc, 3 are candidate transitional 
objects, which is a high rate for a region with such a young age 
(2 Myr), pointing to fast inner disc dispersal. However, due to the 
lack of data between 3.5 and 8 fim, we cannot further discuss the 
differences in inner disc structure; additional observations could 
shed more light on this topic. For the four non-detected mem- 
bers, we derived upper limits, to put some constraints on their 
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Table 6. The logarithm of the IR excess over the photosphere at 
10, 24, 30 and 70 yum. Upper limits are indicated with '<'. 
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Fig. 11. The logarithm of the IR excess at selected wavelengths, 
plotted in function of the effective temperature of the central star. 
Upper limits are indicated with arrows. The warmer the star, the 
larger the infrared excess - apart for the warmest object, MBM 
12-5, with a temperature of 4660 K. 



disc properties. In all 4 cases, the upper limit at 10 micron does 
not allow for an IR excess to be present, suggesting that, if these 
objects do have a disc - what cannot be excluded from the upper 
limits derived at longer wavelengths - these discs must have in- 
ner holes. In Table|6] we list the excess luminosity over the stel- 
lar photosphere; the data are visualised in Fig.Qjj where we plot 
the excesses versus the effective temperature. Those sources that 
have the largest excesses, also have the highest temperatures, al- 
though there is no linear relation between both quantities. 

4.3.2. The influence of companions 

The radiation of cold, outer discs was detected by Hogerheijde 
et al. d2003l > at 450 and 850 /an for 4 stars in MBM 12: LkHa 
262 (MBM 12-2), LkHa 263 ABC (MBM 12-3), LkHa 264 A 
(MBM 12-4) and S18 ABab (MBM 12-12). These authors de- 
rived disc masses between 0.005 and 0.23 xl0~ 2 M© for these 4 
sources, and noticed a relation between the cold disc masses (as 
derived from the submm photometry) and the distance of a com- 
panion: the 2 stars with separations of 100-200 AU have much 
smaller disc masses than the 2 stars with separations of 2000- 
4000 AU. Artymowicz & Lubow (1994) discussed disc trunca- 
tion by companions, and showed that the tidal limit for a disc 
around a member of a binary is ~ 0.4 times the separation. For 
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objects with a not too wide companion, relative to their disc size, 
this means that the companion will effectively truncate their disc. 
Jensen et al. ( 1996) showed with millimetre observations that in- 
deed, objects with separations less than ~ 100 AU have smaller 
disc masses, when compared to their single counterparts, or to 
those objects which have only wider companions. 

We now want to see what the influence of companions is on 
the 'middle' region of the disc, that can be observed in the mid- 
IR. Therefore, we calculated the IR excess at 30 micron for 2 
groups we defined in our sample: 1) the 'binaries', including bi- 
naries with projected separations smaller than 400 AU (MBM 
12-3, 5, 10 and 12), and 2) the 'singles', including single stars 
and binaries with larger projected separations (MBM 12-2, 4 and 
6). This gives average excess values of F30 / F*^o ~ 25 for the bi- 
naries, and f3o/^*,3o ~ 80 for the singles, indicating that the 
discs around the binaries are indeed influenced by their compan- 
ions, as their excess is a factor 3 smaller than that of the single 
stars. 

Bouwman et al. (2006) studied a sample of T Tauri stars 
in the 8 Myr-old cluster 77 Cha, and found that only 1 out of 6 
known or suspected binaries (all with projected separations < 20 
AU) retain a protoplanetary disc. Conversely, 7 out of 9 stars be- 
lieved to be single have discs, suggesting that binary discs have 
shorter lifetimes. Our MBM 12 sample does not include such 
close companions around the disc-bearing stars, but nevertheless 
shows a dependence of its dust characteristics on the presence of 
a companion. It is interesting to see that the two groups (binaries 
and singles) already differ in disc properties at such a young age, 
even when there are no very close companions involved. 

4.3.3. Disc flaring 

The classical model of a T Tauri disc is one in which the op- 
tically thick midplane (containing the bulk of the disc mass) 
is surrounded by an optically thin disc atmosphere, which is 
hot and source of the infrared emission features. The disc can 
flare under hydrostatical equilibrium, and when the gas and dust 
are well-mixed, the turbulent motion of the gas will support the 
dust grains in the disc atmosphere against gravitational settling 
towards the disc midplane. In the disc evolution model from 
Dullemond & Dominik (2004), grains located in the atmosphere 
of a flared disc that grow can no longer be supported by the gas. 
As a consequence, the larger grains will settle towards the disc 
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Fig. 12. Disc properties related to the amorphous grain size, 
derived from the warmer region. Top panel: Degree of flaring 
(derived from the flux ratio at 24 and 8 yum) versus the mass- 
averaged grain size. The larger the grains, the smaller the flar- 
ing, what points to dust settling when grains have grown. Lower 
panel: Relation between the equivalent width of Ha and the 
mass-averaged size of the amorphous grains. The trend suggests 
that, the larger the accretion rate, the larger the grains in the disc 
atmosphere are. 



midplane, hereby reducing the amount of flaring in the disc. In 
Fig. [12] we plot the ratio of the fluxes at 24 and 8 yum versus the 
derived grain size of the amorphous silicates. This flux ratio can 
be used as a proxy for the degree of flaring in a disc (e.g. Apai 
et al. 120051 Scholz et al. |2007| for brown dwarfs and Furlan et al. 
2005 for TTS): the larger the ratio, the more flared the disc. It is 
clear that the more flared the MBM 12 discs are, the smaller the 
derived silicate grain sizes, in agreement with the dust settling 
model from Dullemond & Dominik d20041 > . 

4.3.4. Accretion rate 

We only have two weak-line T Tauri stars in our detected sample: 
MBM 12-1, which has no excess emission, and MBM 12-10, 
which has an IR excess due to warm dust. We thus see accretion 
in all but one of the objects which have an IR excess. Sicilia- 
Aguilar et al. (2007) showed a particular relation between the 
accretion rate derived from the U-band emission, and the grain 
sized derived from the 10 fim feature: they found that objects 
with higher accretion rates have, in general, larger dust grains. 
This suggest that turbulence supports large dust grains against 
settling, and that when the accretion weakens, it will become ap- 
parent in the observed grain size. In Fig.[T2l we plot the equiv- 
alent width of the Ha line in function of the derived grain size, 
and we find a similar tendency for our T Tauri stars - apart from 
object MBM 12-10, which is a WTTS. It is thus tempting to as- 



sume that, the larger turbulence is, the easier it is to support large 
grains against settling towards the disc midplane. Caution should 
be taken, however, given 1) the variability of the Ha line, and 
2) the spectral dependency of the equivalent width for a given 
accretion rate. Furthermore, some weak-line TTS - based on 
the Ha equivalenth width - were found to anyway be accreting, 
through the analysis of the width of their broad Ha line (Sicilia- 
Aguilar et al. 2006). Therefore, additional studies, where the line 
is observed at several epochs and with higher spectral resolution 
(so that the accretion rate can be derived from the 10 % veloc- 
ity width instead of the equivalent width) are required to better 
understand this relation. 



5. Conclusions 

We presented Spitzer spectroscopy and photometry for the com- 
plete sample of T Tauri stars in the 2 Myr-old star forming cloud 
MBM 12. From the 12 objects observed, 8 were detected and 
could be studied more in detail in this paper. We composed spec- 
tral energy distributions for our sample, and found that only 1 
object shows no excess at all, and is probably a disc-less T Tauri 
star. The other objects show excesses that can be attributed to the 
presence of a protoplanetary disc. Of the 7 IR excess sources, 3 
have candidate transitional discs. 

We analysed the properties of the warm dust (visible through 
their solid-state bands in the infrared) with the two-layer tem- 
perature distribution method from Juhasz et al. (2009) and de- 
rived the mass fraction of the different species, and their mass- 
averaged sizes. We showed that the dust components seen in 
other young stars can very well explain the spectral features 
observed in MBM 12 discs. We also found evidence for grain 
growth: the shape and the strength of the 10 fim feature shows 
that grains with submicron sizes, as well as grains of sizes up 
to a few microns are present in the discs in varying amounts. 
The presence of crystalline silicates demonstrate that high- 
temperature processes take place early around other stars, con- 
sistent with meteoritic evidence from the young solar system. In 
a next step, we related the derived properties with the stellar and 
disc properties, and found: 

1. The later the spectral type of the central object, the weaker 
the 10 fim feature, implying larger dust grains. This is most 
likely a luminosity effect, as in less luminous stars, the ra- 
dial distance observed at 10 fim is smaller, hence the density 
larger and, as a consequence, faster grain growth - it does not 
need to imply a real difference in the size distribution 

2. A spatial gradient in grain growth: the mass-averaged grain 
sizes are larger in the warmer region than in the cooler re- 
gion, which also can be related to a difference in density 

3. There is a large spread in the degree of crystallisation (be- 
tween 3 and > 30%), despite a similar age for the objects; 
the crystallisation fraction is independent of spectral type 

4. The crystallisation is also independent of grain growth: the 
mass fraction in crystalline grains nor the size of the crys- 
talline grains are correlated with the size of the amorphous 
silicates. However, the largest crystalline fractions are ob- 
served for those sources with largest grain sizes 

5. A spatial gradient in crystalline composition, as observed 
in the forsterite to enstatite mass ratio: forsterite is much 
more dominating in the outer than in the inner disc re- 
gion. Enstatite is more abundant in those sources with larger 
grains, and in regions where the density is higher. This sug- 
gests that its formation process needs higher densities than 
forsterite 
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6. The presence of a spatial gradient in both grain sizes and 
composition indicates that radial mixing is not a very effi- 
cient process in these protoplanetary discs 

7. Companions, even at distances larger than 70 AU, influence 
disc evolution: the IR excess of single sources is a factor 3 
higher than that of binaries. 

8. Disc flaring and grain size are mildly related: the more the 
discs flare, the smaller the grains in their atmosphere 

9. There is tendency for objects with a stronger Ha line to have 
larger grains, suggesting that more heavily accreting (more 
turbulent) sources can support larger grains in their disc at- 
mosphere more efficiently 

We are aware that these findings are based on a small sample, 
and that it should be confronted with a larger database of well- 
known T Tauri stars, for which infrared spectra of high qual- 
ity are available. This study is part of a larger programme with 
Spitzer on discs in young clusters. The low-mass M-type stars 
in the Coronet cluster have already been discussed in Sicilia- 
Aguilar et al. (2008 ), while T Tauri stars in the r\ Cha cluster 
will be presented by Sicilia-Aguilar et al. (in preparation) and T 
Tauri stars in the e Cha cluster by Fang et al. (in preparation). 
From a comparison of the results in the different clusters, we 
will be able to better determine the relations between the differ- 
ent stellar and dust properties, so that they can be incorporated 
in theoretical (evolutionary) disc models of young objects. 
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Table 7. Parameters involved in the TLTD model 



Parameter 


Meaning 


Fy 


Observed flux 




Total continuum flux 


d 


Distance to the source 


R, 


Stellar radius 


T, 


Stellar temperature 


D 


Contribution of the star to the total flux 


Di 


Contribution of the inner rim to the total flux 


D 2 


Contribution of the midplane to the total flux 


Di,j 


Mass contribution of species i with size j 


T a , out 


Lowest temperature in the disk atmosphere 


Ta, in 


Highest temperature in the disk atmosphere 


J* out 


Lowest temperature in the inner rim 


T 


Highest temperature in the inner rim 


T 

1 tn, out 


Lowest temperature in the disk midplane 


Tm, in 


Highest temperature in the disk midplane 




Mass absorption coefficient of species i with size j 


qr 


Power exponent of the temperature distribution (as a function of radius) in the rim 


qa 


Power exponent of the temperature distribution (as a function of radius) in the disk atmosphere 


qm 


Power exponent of the temperature distribution (as a function of radius) in the disk midplane 
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Table 8. For each object, we list the reduced^ 2 of the fit, and the mass fraction of each species for the 3 different sizes, used to fit 
the shorter and longer wavelengths regions (upper and lower part of the table, respectively). We only show mass fractions that are 
higher than 5%. 



MBM 12-2 MBM 12-3 MBM 12-4 MBM 12-5 MBM 12-6 MBM 12-10 MBM 12-12 



7-17 micron range 


x 2 


12.2 


5.7 


9.0 


9.7 


4.9 


3.4 


5.6 




Size 






Mass fractions (in percentage) 






Amorphous silicate 


0.1 yum 


13.2 C&J) 


- 


48.2 C£?) 


43.4 O 


- 


- 


- 


(Olivine Type) 


1.5 fim 


- 


- 








- 


- 




6.0 fim 


- 


29.7 es£) 






8.8 Cf") 


17.6 (^) 


53.3 £11) 


Amorphous silicate 


0.1 yum 






15.9 ( +3 2 l) 










(Pyroxene Type) 


1.5 fim 


20.2 (*>■») 


35.2 ( + _ 9 7 D 


16.8 eg) 


24.0 O 




13.8 CH) 






6.0 yum 


55-7 O 


19.8 C1U) 




29.0 C£) 


32.0 


62-0 O 


10.9 CJJ) 


.TUIalCIllC 


ft 1 nm 


















1.5 fim 


















6.0 fim 


_ 


_ 








_ 


_ 


Enstatite 


0.1 


_ 


_ 








_ 


_ 




1.5 fim 


- 


— 








- 


- 




6.0 fim 


9.0 O 


_ 


10.6 cii) 




40.1 C|^) 


_ 


29.1 ci# 


Silica 


0.1 fim 




_ 








_ 






1.5 fim 


- 


— 








— 


— 




6.0 fim 


_ 


7.0 Cg) 






5.4 C««) 


_ 


_ 


17-37 micron range 


x 2 


12.3 


8.3 


14.0 


14.2 


5.9 


3.6 


5.7 


kj UL'LHyJ 


Size 






Mass fractions (in percentage) 






Amorphous silicate 


0.1 fim 


65.2 (™) 


_ 


51.1 qp 








30.5 C™) 


(Olivine Type) 


1.5 fim 




68.8 CJJ) 


6.4 (f 3 4 ) 


36.7 






32.7 Clg) 




6.0 fim 


- 














Amorphous silicate 


0.1 fim 


30.4 C_lf) 


21.0 C||) 


33.0 


47.6(^ 8 6 ) 


8.1 Cf, 2 ) 


41.7 Cul) 


31. 5 C£f) 


(Pyroxene Type) 


1.5 yum 










79.4 (!*») 








6.0 yum 
















Forsterite 


0.1 fim 










5.8 O 








1.5 yum 


















6.0 yum 
















Enstatite 


0.1 yum 


















1.5 yum 


















6.0 yum 
















Silica 


0.1 yum 


















1.5 yum 


















6.0 yum 






6.4 O 


7.1 O 




8.7 CtD 





